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Abstract — A novel general ready-to-use bit-error rate (BER) 
expression for one-dimensional constellations is developed. The 
BER analysis is performed for bit patterns that form a labeling. 
The number of patterns for equally spaced A/-PAM constellations 
with different BER is analyzed. 

I. Introduction and Motivation 

Current wireless communication systems are based on the 
bit-interleaved coded modulation (BICM) paradigm introduced 
in Qj and later studied in (O, 0. One key element in these 
systems is the calculation of logarithmic likelihood ratios 
(LLR, also known as L-values) for the received bits, which are 
passed to the channel decoder. The coded performance analysis 
of such systems is generally not straightforward, and is usually 
carried out either numerically by Monte-Carlo simulation, or 
in terms of lower and upper bounds (21 Sec. 4], (5] Ch. 4]. 
The calculation of LLRs is crucial also in many other coded 
systems. In this paper, we analyze the uncoded performance 
over the additive white Gaussian noise (AWGN) channel. 

A symbol-based demodulator (SD) is the most natural way 
of decoding symbols transmitted through the channel. This 
approach is optimal in terms of symbol-error rate (SER). The 
bit-error rate (BER) performance of the SD is well documented 
in literature, e.g. [4, Ch. 5], Gl Ch. 10], 0-lHT] and references 
therein. On the other hand, in a coded system, such as BICM, 
soft or hard information on the received bits is passed to 
the decoder, and thus, bit-wise decisions are more relevant 
than symbol-wise decisions. The optimal bit-wise demodulator 
(BD) minimizing the BER implies the calculation of (exact) L- 
values for the received bits. The uncoded performance of such 
a demodulator has been studied in fl2l . where closed-form 
expressions for the BER for 4-ary pulse amplitude modulation 
(PAM) with the binary reflected Gray code (BRGC) l9l. |[T3l 
|[T4l are presented. Due to the complexity of the BD, the 
calculation of L-values in practical systems is usually done 
based on the so-called max-log approximation lfl5l eq. (5)], 
lfl6l eq. (1)]. We call this demodulator the approximate BD 
(ABD). The three above demodulators (SD, BD, and ABD) 
have been recently numerically compared from a mutual 
information point of view in ifTTl for multiple-input multiple- 
output BICM systems. 
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In this paper, we prove the equivalence of the SD and 
the ABD in terms of uncoded BER for any constellation 
and labeling. Due to this equivalency, we go on and study 
the ABD for one-dimensional constellations. To this end, 
we introduce a novel ready-to-use BER expression valid for 
any one-dimensional constellation and binary labeling. The 
analysis is performed for bit-patterns that form a labeling. 

II. Preliminaries 

A. Notation Convention 

The following notation is used throughout the paper. Lower- 
case letters x denote real scalars and boldface letters x denote 
a row vector of scalars. Blackboard bold letters X denote 
matrices with elements Xij in the ith row and the jth column 
and (-) T denotes transposition. Calligraphic capital letters X 
denote sets, where the set of real numbers is denoted by 1Z. 
The binary complement of x 6 {0, 1} is denoted by x — 1—x. 
Binary addition (exclusive-OR) of two bits a and b is denoted 
by a © b. Random variables are denoted by capital letters X 
and probabilities by Pr{-}. The Gaussian Q-function is defined 

as Q(x) = -4= L°°e--T dt. 



B. System Model 

In this paper we analyze a system where a vector of 
binary data b = [bi , . . . , b m ] is fed to a modulator. The 
modulator carries out a one-to-one mapping from b to one 
of the M constellation points x G X = {si, . . . , sm}, where 
s i < s 2 < • ■ ■ < Sm> f° r transmission over the physical 
channel, where M — 2 m . The modulator is defined as the 
function $ : {0,l} m — > X with a corresponding inverse 
function : A" — > {0, l} m . 

For PAM constellations, Sj = -d(Al — 2i + 1),£ = 
l,...,M, where d = y/i/JW — 1) to normalize the constel- 
lation to unit average energy, i.e., E s = jj Yli=i s 1 — 1- We 
assume the bits to be independent and identically distributed 
(i.i.d.) with Pi{B 3 = u} = 0.5, Vj and u G {0, 1}, and thus, 
the symbols are equiprobable, i.e., Pr{X = s;} = 1/M, Mi. 

The modulator is defined by the constellation and its binary 
labeling. A binary labeling is specified by the matrix C = 
[c[, . . . , cjj] 7 of dimensions M by m, where the ith row 
Ci = [ci,ii ■ ■ ■ i c i,m\ is the binary label of the constellation 
point Si, i.e., $(ci) = S{. 



In this paper we consider a discrete time memoryless 
AWGN channel with output y = x + rj, where x G X and the 
noise sample T) is a zero-mean Gaussian random variable with 
variance No/2. The conditional PDF of the channel output is 
given by 

PY \x{y\x) = ^- p[v - x) \ (1) 

where the average signal to noise ratio (SNR) is defined as 
p 4 E s /N = 1/iVo. 

The observation y is used by the demodulator to decide 
on the transmitted binary sequence, i.e., to produce b = 
[bi, ■ ■ -,b m ]. 

C. Demodulators 

The SD makes a hard decision on the transmitted symbol 
and returns the length-m binary label of that symbol, i.e., 
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6°" 4 $ 1 ( argmin (y - x) 2 

xex 



(2) 



The SD in (O is optimal in terms of minimizing the SER, but 
it does not necessarily minimize the BER. 

To minimize the BER the optimal BD should be used. The 
BD calculates (a posteriori) L-values for the m bits based on 
the observation y, i.e., 

Pi{Bj =l\Y = y} 



hiy) = 1 °? 



log 



Vv{Bj =0\Y = y} 



YxeXj.o ' 



-p(y-x) 2 



(3) 



(4) 



where j = 1, . . . , m and Xj jU = {si G X : Cfj = u, Vi}. To 
pass from (O to d4j Bayes' rule was used together with the 
i.i.d. assumption of the bits and the conditional PDF in ([1). 

The implementation of the BD in its exact form © is 
complicated, especially for large constellations, as it requires 
calculation of the logarithm of a sum of exponentials. To 
overcome this problem, approximations are usually used in 
practice. The most common approximation is the so-called 
max-log approximation (logJ2 i e Ai w max^ Ai) 0] eq- (3.2)], 
E eq. (9)], lEl eq. (5)], Ql eq. (8)], which used in © gives 



l Av) = p 



min (y — x) 2 — min (y — x) 2 

x&Xjfi x£Xji 



(5) 



The ABD is defined as the demodulator that applies the 
following decision rule 



? ABD 



1 if/,(y)>0, 
otherwise. 



(6) 



The next theorem gives proof for the equivalence of the 
SD and the ABD. This was mentioned in flj] Sec. IV-A], 
however, no proof was given there. 

Theorem 1: For any p, X, and C, bj = b^ BB for all j — 
1, . . . ,m. 

Proof: Combining (0 and ([5]), the decision rule for the 



ABD can be written as 

£abd = I 1 - mirix£*3,o (y - x) 2 > mm xl£Xj l (y ~ x) 2 , 
3 [0, min xeX] Q (y - x) 2 < iam xeXj l (y - x) 2 , 



which can be simplified to 



lABD 



argmin < min {y — xY 

nG{0,l} I xeX i 



(7) 



Since min ue{0il } {min x£Xj u (y - x) 2 } = min^* (y - x) 2 
for any X, p, and C, the symbol found by the ABD in (0 
will always be the closest x G X to y in terms of Euclidean 
distance (ED), regardless of the bit position j. This is the same 
rule used in (|2]), which completes the proof. □ 
Theorem \T\ states that the SD and the ABD are equivalent 
and optimal in terms of minimizing the SER for any constella- 
tiorfl and any labeling. Because of this, from now on we only 
consider the ABD. 

III. BER for One-Dimensional Constellations 
The BER for a given labeling C can be expressed as 



Pi 



1 m 
m ^-^ J 



where the BER for the jth bit position Pj 
bj\Bj = bj} can be written as 



(8) 



1 M 

i=l 



Si} 



(9) 



using the law of total probability. The BER for the jth bit 
position Pj depends only on the subconstellations Xj and 
Xj i (cf. ©-(fSJl), i.e., on the jth column of C, such that Pj = 

P([c hj ,...,C M ,j})- 

We define a bit pattern (or simply pattern) as a length-M 
binary vector p = [pi, . . . ,Pm] G {0, 1} m with Hamming 
weight M/2. The labeling C can now be defined by m 
patterns, each corresponding to one column of C. We index 
the patterns as p w with w bein^ the decimal representation 
of the vector p, i.e., w = Xa=i 2 M ~ l pi. For example, for 
M = 4, the pattern [0, 1, 0, 1] is indexed as p 5 (cf. Table [Q). 
The BER for the labeling C does not depend on the order 
of its columns, and thus, the BER for the labeling C is fully 
determined by a set of m patterns W = {wi, . . . , w m }. 

Based on the previous discussion, from now on we concen- 
trate our analysis only on patterns (and not on labelings), i.e., 
on the function P(p), however, to simplify the notation, the 
dependency on the pattern will be omitted. 

To analyze the BER of a pattern (PBER), the observation 
space 1Z is split into two disjoint decision regions, i.e., To = 
{y G K : b = 0} and Ti = {y G K : b = 1} such that 
r U Ti = K. 



'The proof of TheoremfT]was given for one-dimensional constellations only, 
however, its extension to any multi-dimensional constellation is straightfor- 
ward. 



Using the definition of Tq and Ti, the PBER for the pattern 
p can be rewritten as 



1 



M 



i=l 



(10) 



By expressing P as in ( fTOl l, it is clear that the PBER in (0 
can be calculated using the decision regions r and Ti only, 
as opposed to alternative approaches where ( fTOt is expressed 
in terms of the PDF of the L-values (cf. [19 eq. (19)], (20] 
Sec. IV]). 

Decision thresholds (or simply thresholds), denoted by j3 k , 
where k — 1,2,... stands for the index of the threshold, 
are defined as the points that separate the decision regions 
for zeros and ones, and thus, they fully determine the PBER 
in ([Tol l. The thresholds for the ABD are the midpoints between 
the constellation points labeled with different bits, which 
follows directly from (0. 

The BER expression for the ABD and an M-PAM con- 
stellation with any labeling is well known and can be found 
in lfT4l eq. (21)]. The PBER expression can easily be obtained 
in a similar way. In the following theorem, we generalize the 
result in lfl4l eq. (21)] to non-equally spaced constellations and 
derive a general PBER expression for any one-dimensional 
constellation. 

Theorem 2: The PBER for the ABD using an arbitrary one- 
dimensional constellation with a pattern p can be expressed as 



M M-l 



p =2 + MT,T,^kCl((f3k- Sl )V^p)^ (ID 

i=l k=l 

where /3 fc = Sfc+ ^ +1 , k= 1, . . . , M-l and g i>k G {0, ±1} is 



9i,k = (Pk+i ^ 2pi). 



(12) 



Proof: Let Vi jk be the following conditional probabilities 

v i ,i=Pt{Y<0 1 \X = a i } 

= 1 - Q {{fix - Si) ^27) , (13) 
v i , k ='PT{fik-i<Y<0 k \X = a i } 

= Q (09*_ x - Si)y/2p) - Q ((& - 8i)y/Tp) , (14) 
Vi,M = Pt{Pm-i < Y\X = Si } 

= Q ((Pm-i — Si)\/%?J , (15) 

where i = 1, . . . , M, k = 2, . . . ,M - 1, and /3 k = Sfc+ 2 Sfc+1 
for k = 1, . . . , M — 1. The PBER in (fTOb can be rewritten as 

1 M 



M M 



i=i k=i 

where e i>k = Pi®Pk e {0, 1}. 



(16) 



Using ffT3T> — (fT3T> the PBER in (fTol l can be expressed as 

M M M 



P =M 



i=l k=2 
M M-l 



i=i k =i 

M M-l 



2 + M E E ^ i,fe+1 ~ &i ' k ^ ~ s ^^p) ' 

(17) 



i=i k =i 



where 2~2iLi e «.i — 2~2t=i Pi ®Pi = M/2 was used. To obtain 
the expression in ( fTTT i. we express e^fe+i — e^fc in fill as 

e ijk+1 - e ijk = p k+ \ 8 Pi - p k ® Pi (18) 
= (pfc+i-j»*)(l-2Pi). (19) 

where the identity p L © pj = pipj + piPj was used together 
with pi = 1 - p t . 

The threshold /3 k between the constellation points labeled 
with the same bit does not affect the PBER in ( fTTT ) as g ik — 
0, Vi in <[T2). □ 

Remark 1: Theorem |2] gives an expression for the PBER 
for the ABD. However, (fTTT i can be used for calculating the 
PBER when the thresholds (5 k are not midpoints or, moreover, 
when they are dependent on the SNR, for example, when the 
BD is used. Analytical expressions for thresholds for the BD 
are in general unknown. 

To illustrate Remark Q] consider 8-PAM labeled by 
the BRGC, which is formed by the patterns p 15 = 
[0,0,0,0,1,1,1,1], p 60 = [0,0,1,1,1,1,0,0], and p W2 = 
[0,1,1,0,0,1,1,0]. From <[TTJ-jT2j, whenever g lM = 0, the 
value of f3 k does not influence the PBER and can be set to 
any value. The thresholds for gi, k ^ can be numerically cal- 
culated setting the L-value in (@]i to zero. The obtained results 
are shown in Fig. [T] Using these thresholds in (fTTT i and ©, 
the BER for the patterns and for the BRGC are calculated. 
The results for the BD and the ABD are presented in Fig. [2] 
and show no notable difference between the demodulators for 
p > dB. 

IV. BER FOR M-PAM 

In this section, we study the BER for equally spaced M- 
PAM constellations. We concentrate on classifying patterns 
and comparing their performance. For M-PAM, ( fTTT i can be 
expressed as a bit-wise version of [14, eq. (21)]: 



M-l 



(20) 



where 



M-l 



hi = 22 ( Pk+1 _ Pk ^ 1 ~ 2 Pfe+l-n) 

- (pk+2-n - Pfc+1-„)(1 - 2p k+ i). (21) 



k=n 
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Fig. 1. Thresholds for 8-PAM with different patterns vs. SNR. Due to the 
symmetry of the patterns the thresholds are symmetric with respect to zero. 
Only positive thresholds are shown. Squares represent the constellation points. 
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Fig. 2. The BER for 8-PAM with patterns p 15 , p 60 , p 102 , and the BRGC. 
Solid lines correspond to the BD and dashed lines correspond to the ABD. 



One direct consequence of ( f20~b is that the vector a = 
[ai, . . . , <Xm-i] with a n given by d2lT i completely defines 
the performance of the ABD for ilf-PAM and allows us to 
compare the performance of different patterns. From J20b . 
the PBER for high SNR can be predicted by the coefficient 
multiplying the Q-function with the smallest argument, that is, 
a\. If for two patterns the coefficients are identical, the next 
coefficients 0,2 are checked, and so on. 

We observe that, for instance, for 4-PAM, the pattern 



p 5 = [0,1,0,1] and the pattern p 10 = [1,0,1,0] have 
identical PBER performance because of the symmetry of the 
constellation. It is therefore interesting to find all the patterns 
with different performance. This will allow us to predict the 
performance of any possible labeling. We therefore group all 
the patterns with identical performance into one class. The 
next theorem gives a closed form expression for the number 
of classes for length- M patterns. 

Theorem 3: For M-PAM, all the length-M patterns can be 
grouped into Q classes, where the patterns within each class 
have identical PBER, and 

Q = i(,M/ 2 ) + 0+ 2M/2 )- (22) 

Proof: We define two operations that can be applied to 
a pattern that will be used in the proof. A reflection of p is 
defined as p' = refl(p) with p[ L — pM+i-i for i = 1, . . . , M. 
An inversion of p is defined as p' = inv(p) with p\ = pi 
for i — 1, . . . , M. Both these functions are self-inverse, i.e., 
p = refl(refl(p)) and p — inv(inv(p)), and they commute, 
i.e., refl(inv(p)) = inv(refl(p)). Note also that for any pattern 
p, we have that p ^ inv(p). 

We introduce three special types of patterns. The pattern 
p is said to be reflected (RE) if refl(p) = p, the pattern p 
is said to be anti- reflected (ARE) if inv(refl(p)) = p, and 
the pattern p is called asymmetric (ASY) if it is neither RE 
nor ARE. For example, p 60 — [0, 0, 1, 1, 1, 1, 0, 0] is an RE 
pattern, p 43 = [0,0,1,0,1,0,1,1] is an ARE pattern, and 
p 216 = [1,1,0,1,1,0,0,0] is an ASY pattern. 

From d9t-([T0ii, we note that the PBER is not affected by 
reflections and/or inversion of the patterns, since the PBER 
is averaged over both transmitted zeros and ones. Because of 
this, we group all patterns that are connected via reflection 
or inversion into one class of patterns that all have identical 
PBER. Each class contains either two patterns (p and inv(p) 
because p ^ inv(p), Vp) or four patterns (p, inv(p), refl(p), 
and inv(rcfl(p))). 

Any pattern p must contain M/2 zeros and M/2 ones, 
hence, the total number of patterns is equal to ( A *y 2 )- For 
a pattern to be RE, pi = pM-i+i, i.e., the positions of the 
M/4 ones in [pi, ■ ■ ■ ,Pm/2] fully describe the pattern, and 
thus, the number of RE patterns is (jj^). There are two 
members in every class of RE patterns, p = rcfl(p) and 
inv(refl(p)) = inv(p), which gives 5(^/4) classes. 

For a pattern to be ARE, pi = pM-i+i, i.e., the positions of 
the ones in [pi, . . . ,Pm/2] fully describe the pattern, where the 
number of ones in [pi, . . . ,Pm /2I is between and M/2. From 
that, it follows that there are 2 M / 2 ARE patterns. There are 
two members in every class of ARE patterns (p = inv(rcfl(p)) 
and rcfl(p) = inv(p)), which gives 2 A/ / 2_1 classes. 

All the remaining classes include only ASY patterns. The 
number of ASY patterns can be obtained by subtracting the 
number of RE and ARE patterns from the total number of 
patterns. There are four patterns in each class, as p ^ rcfl(p) 
and p ^ refL(inv(p)) (or equivalently, rcfl(p) 7^ inv(p)). 
Using this, the total number of classes in ( f22l ) is obtained 
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(b) 16-PAM 

Fig. 3. The PBER for the patterns for 8-PAM and 16-PAM. All the curves 
merge into M — 1 groups at high SNR as predicted by Remark [2] 



as sum of classes of RE, ARE, and ASY patterns. □ 
For example, Theorem [3] states that there are 3 classes of 
patterns for 4-PAM, 23 classes for 8-PAM, and 3299 classes 
for 16-PAM. The PBER for 8-PAM and 16-PAM for all the 
patterns is shown in Fig. [3] All the classes of patterns for 
4-PAM and 8-PAM are shown in the first and the second 
parts of Table U respectively. For each class, Table [Q shows 
the representative of the class p, the decimal indices of class 
members w (the index of the representative is shown with 
boldface), and the vector a that defines the PBER. The patterns 
are ordered from best to worst PBER at high SNR. 

Remark 2: The element a\ in d2"TT i is equal to twice the 



TABLE I 

Classes of patterns for 4-PAM and 8-PAM with their 

CORRESPONDING REPRESENTATIVES p , DECIMAL REPRESENTATIONS OF 
THE PATTERNS W, AND VECTORS a DEFINING THEIR PBER 



p 


w 


a 


[0,0, 1, 1J 


3 12 


[2, 2, 0] 


[0,1,1,0] 


6 9 


[4, 2, -2] 


[0,1,0,1] 


5 10 


[6,-4, 2] 




0, 0, 0, 0, 1, 1, 1, 1 


ID Z4U 




2, 2, 2, 2, 0, 0, 






0, 0, 0, 1, 1, 1, 1, 


Qn i on 1 1 x. oof; 

oU 1ZU loD ZZD 




4, 3, 3, 2, —2, —1, — 1 






0, 0, 1, 1, 1, 1, 0, 


fin i 




4, 4, 2, 2, —2, —2, 






0, 0, 0, 1, 0, 1, 1, 1 


ZO zoz 




6, —2, 2, 0, 2, 0, 






0, 0, 0, 1, 1, 1, 0, 1 


£i\3 I 1 lo4 ZZO 




6, 1, 2, —3, 1, 0, 1 






0, 0, 0, 1, 1, 0, 1, 1 


97" 91 R 998 
z i oy ZIO zZo 




6, 2, —3, 1, 1, 1, 






0, 1, 1, 1, 0, 0, 0, 1 


1 1 O 1 A ry 

J.J.O 14Z 




6, 4, 4, —4, -2, -2, 2 






0, 0, 1, 1, 1, 0, 0, 1 


o / yy loo iyo 




6, 5, 0, -3, -3, 2, 1 






U, U, 1, 1, U, U, 1, 1 


E -| OCX A 
O.L ZU4 




0, 0, —4, —4, 2, 2, U 






0,0,1,0,1,1,1,0 


4to no ioy zuy 




8,-1,2,-1, 3,-2,-1 






0,0,1,1,1,0,1,0] 


09 1 fi'^ 1 07 

do yz ido ly i 




8, -1,3, -2, 2, -1, -1 






0,1,0,0,1,1,1,0 


7Q 11/1 1/11 1 11 
I o 1 14 141 111 




8, 2, -1, -1, -1,3, -2 






0,0,1,1,0,1,1,0 


t^d 1 08 1 A7 901 




8, 3,-6, 3, 3,-2,-1 






0,1,1,0,0,1,1,0 


102 153 




8, 6,-6,-4, 4, 2,-2 






0,0,1,0,1,0,1,1] 


43 212 




10,-6, 4,-2, 0, 2, 






0,0,1,0,1,1,0, 1 


45 75 180 210 




10,-3,-3, 6,-4, 1, 1 






0,0,1,1,0,1,0, 1 


53 83 172 202 




10,-3, 1, 0,-2, 1, 1 






0,1,0,0,1,1,0, 1 


77 178 




10, 0,-6, 2, 4,-4, 2 






0,1,1,0,1,0,0, 1] 


105 150 




10, 0,-4, 6,-4,-2, 2 






0,1,0,1,1,0,0, 1J 


89 101 154 166 




10, 0,-3, 1, 1,-3, 2 






0,1,0,1,1,0,1,0] 


90 165 




12,-6, 0, 6,-6, 4,-2 






0,1,0,1,0,1,1,0] 


86 106 149 169 




12, -6,3, -1, -1, 3, -2 






0,1,0,1,0,1,0, 1] 


85 170 




14,-12, 10, -8,6,-4,2 





number of pairs of constellation points at minimum ED whose 
bits are different (for a given pattern). Using this, it can be 
shown that for il/-PAM there are M — 1 different values of 
a%. This means that the PBER of all the patterns merge into 
M — 1 groups at high SNR. For example, for 8-PAM and 16- 
PAM the number of groups of patterns at high SNR is 7 and 
15, respectively, as illustrated by Fig. [3] 

Using © and (O, the average BER for M-PAM with 
labeling C can be expressed as lfl4l eq. (21)]: 

M-l 

p c = — ; 77 E «„Q ((2n - l)d^2~ P ) , (23) 

n—l 

where a = [a>i, . . . , aju-i] is the sum of vectors a for the 
m patterns used in C. The value of a n is a scaled version 
of the so-called differential average distance spectrum 5(n, A) 
in HI eq. 21], i.e., a n = 2M8(n, A). 

Remark 3: The value of ct\ corresponds to twice the sum of 
Hamming distances between binary labelings of constellation 
points at minimum ED. It can be shown that = 2m(M — 
1) — a\, where A^, was recently shown to determine the BICM 
mutual information in the high SNR regime |T2D . 

By listing the vectors a for all the possible labelings for 
8-PAM, we found 12 different ay, which is in agreement with 
the 12 classes of labelings (with different A$) shown in ll2"Tl 
Fig. 2(b)]. The BER for all the labelings for 8-PAM is shown 
in Fig. |U where the 12 classes are visible for high SNR. 

To conclude, we present the vectors a. for 4-PAM and 
8-PAM with some common labelings, including the BRGC, 
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Fig. 4. The BER for all the 460 labelings with different BER for 8-PAM. 



TABLE II 

Some common labelings for 4-PAM and 8-PAM with their 

CORRESPONDING PATTERN INDICES W AND VECTORS a DEFINING THEIR 

BER 



M 


Labeling 


W 


a 


4 


BRGC 


{3,6} 


[6,4,-2] 


4 


NBC 


{3,5} 


[8, -2, 2] 


4 


AG 


{5,6} 


[10,-2,0] 


8 


BRGC 


{15,60, 102} 


[14,12,-2,0,2,0,-2] 


8 


FBC 


{15,60,90} 


[18,0,4, 10, -8,2, -2] 


8 


NBC 


{15,51,85} 


[22,-4,8,-10,8,-2,2] 


8 


BSGC 


{105,60, 102} 


[22,10, -8,4,-2,-2,0] 


8 


AG 


{90, 105,85} 


[36,-18,6,4, -4,-2,2] 



the natural binary labeling (NBC) EH Sec. II-B], the folded 
binary code (FBC) [ 1 1 [22, Sec. II-B], the binary semi-Gray 
code (BSGC) (22l Sec. II-B], and the so-called anti-Gray (AG) 
labeling |23l . These labelings are shown in Table Ull together 
their pattern indices W and vectors a, in the first part for 
4-PAM, and in the second part for 8-PAM. The labelings are 
also ordered from best to worst BER at high SNR. By listing 
the vectors a. for all the possible labelings, we found three 
labelings with different BER for 4-PAM listed in Table HI] For 
8-PAM we found 460 labelings as shown in Fig. |4] 

V. Conclusions 

A novel general expression for the uncoded BER of one- 
dimensional constellations has been introduced. For equally 
spaced A/-PAM constellations, a classification of the patterns 
has been performed and a closed form expression on the 
number of patterns that give different BER has been derived. 
The rule for combining patterns into a labeling remains for 
future investigation. Establishing this rule will allow us to 
define a number of labelings with different BER for an 



arbitrary constellation size M and also a number of groups 
of labelings at high SNR. 

References 

[1] E. Zehavi, "8-PSK trellis codes for a Rayleigh channel," IEEE Trans. 
Commun., vol. 40, no. 3, pp. 927-946. May 1992. 

[2] G. Caire, G. Taricco, and E. Biglieri, "Bit-interleaved coded modula- 
tion," IEEE Trans. Inf. Theory, vol. 44, no. 3, pp. 927-946, May 1998. 

[3] A. Guillen i Fabregas, A. Martinez, and G. Caire, "Bit-interleaved 
coded modulation," Foundations and Trends in Communications and 
Information Theory, vol. 5, no. 1-2, pp. 1-153, 2008. 

[4] J. G. Proakis, Digital Communications, 4th ed. McGraw-Hill, 2000. 

[5] M. K. Simon, S. M. Hinedi, and W. C. Lindsey, Digital Communication 
Techniques: Signed Design and Detection. Prentice Hall, 1995. 

[6] K. Cho and D. Yoon, "On the general BER expression of one- and two- 
dimensional amplitude modulations," IEEE Trans. Commun., vol. 50, 
no. 7, pp. 1074-1080, Jul. 2002. 

[7] P. J. Lee, "Computation of the bit error rate of coherent Af-ary PSK 
with Gray code bit mapping," IEEE Trans. Commun., vol. COM-34, 
no. 5, pp. 488-491, May 1986. 

[8] J. Lassing, E. G. Strom, E. Agrell, and T. Ottosson, "Computation of the 
exact bit-error rate of coherent M -ary PSK with Gray code bit mapping," 
IEEE Trans. Commun., vol. 51, no. 11, pp. 1758-1760, Nov. 2003. 

[9] E. Agrell, J. Lassing, E. G. Strom, and T. Ottosson, "On the optimality 
of the binary reflected Gray code," IEEE Trans. Inf. Theory, vol. 50, 
no. 12, pp. 3170-3182, Dec. 2004. 
[10] J. Lassing, E. G. Strom, E. Agrell, and T. Ottosson, "Unequal bit- 
error protection in coherent M -ary PSK," in IEEE Vehicular Technology 
Conference (VTC-Fall), Orlando, FL, Oct. 2003. 
[11] L. Szczecinski, C. Gonzalez, and S. Aissa, "Exact expression for the 
BER of rectangular QAM with arbitrary constellation mapping," IEEE 
Trans. Commun., vol. 54, no. 3, pp. 389-392, Mar. 2006. 
[12] M. K. Simon and R. Annavajjala, "On the optimality of bit detection of 
certain digital modulations," IEEE Trans. Commun., vol. 53, no. 2, pp. 
299-307, Feb. 1988. 
[13] F. Gray, "Pulse code communications," U. S. Patent 2632058, Mar. 
1953. 

[14] E. Agrell, J. Lassing, E. G. Strom, and T. Ottosson, "Gray coding for 

multilevel constellations in Gaussian noise," IEEE Trans. Inf. Theory, 

vol. 53, no. 1, pp. 224-235, Jan. 2007. 
[15] A. J. Viterbi, "An intuitive justification and a simplified implementation 

of the MAP decoder for convolutional codes," IEEE J. Set Areas 

Commun., vol. 16, no. 2, pp. 260-264, Feb. 1998. 
[16] Ericsson, Motorola, and Nokia, "Link evaluation methods for high 

speed downlink packet access (HSDPA)," TSG-RAN Working Group 

1 Meeting #15, TSGR1#15(00)1093, Tech. Rep., Aug. 2000. 
[17] P. Fertl, J. Jalden, and G. Matz, "Performance assessment of MIMO- 

BICM demodulators based on mutual information," IEEE Trans. Sig. 

Proc, vol. 60, no. 3, pp. 1366-1382, Mar. 2012. 
[18] P. Robertson, E. Villebrun, and P. Hoeher, "A comparison of optimal and 

sub-optimal MAP decoding algorithms operating in the log domain," in 

IEEE International Conference on Communications (ICC), June 1995. 
[19] M. Benjillali, L. Szczecinski, S. Aissa, and C. Gonzalez, "Evaluation of 

bit error rate for packet combining with constellation rearrangement," 

Wiley Journal Wireless Comm. and Mob. Comput., pp. 831-844, Sep. 

2008. 

[20] A. Alvarado, L. Szczecinski, R. Feick, and L. Ahumada, "Distribution of 
L-values in Gray-mapped M 2 -QAM: Closed-form approximations and 
applications," IEEE Trans. Commun., vol. 57, no. 7, pp. 2071-2079, 
July 2009. 

[21] A. Alvarado, F. Brannstrom, and E. Agrell, "High SNR bounds for the 
BICM capacity," in IEEE Information Theory Workshop (ITW), Paraty, 
Brasil, Oct. 2011. 

[22] E. Agrell and A. Alvarado, "Optimal alphabets and binary labelings 
for BICM at low SNR," IEEE Trans. Inf. Theory, vol. 57, no. 10, pp. 
6650-6672, Oct. 2011. 

[23] S. ten Brink, J. Speidel, and R.-H. Yan, "Iterative demapping and de- 
coding for multilevel modulation," in IEEE Global Telecommunications 
Conference (GLOBECOM), Sydney, Australia, Nov. 1998. 



